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ABSTRACT 

Quasi-static speckles are a current limitation to faint companion imaging of bright 
stars. Here we show through simulation and theory that an adaptive pupil mask can 
be used to reduce these speckles and increase the visibility of faint companions. This 
is achieved by placing an adaptive mask in the conjugate pupil plane of the telescope. 
The mask consists of a number of independently controllable elements which can either 
allow the light in the subaperture to pass or block it. This actively changes the shape 
of the telescope pupil and hence the diffraction pattern in the focal plane. By randomly 
blocking subapertures we force the quasi-static speckles to become dynamic. The long 
exposure PSF is then smooth, absent of quasi-static speckles. However, as the PSF 
will now contain a larger halo due to the blocking, the signal to noise ratio (SNR) 
is reduced requiring longer exposure times to detect the companion. For example, in 
the specific case of a faint companion at 5X/D the exposure time to achieve the same 
SNR will be increased by a factor of 1.35. In addition, we show that the visibility 
of companions can be greatly enhanced in comparison to long-exposures, when the 
dark speckle method is applied to short exposure images taken with the adaptive 
pupil mask. We show that the contrast ratio between PSF peak and the halo is then 
increased by a factor of approximately 100 (5 magnitudes), and we detect companions 
11 magnitudes fainter than the star at 5X/D and up to 18 magnitudes fainter at 
22.5X/D. 
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1 INTRODUCTION 

Detecting the faint reflected or self-luminous signal from ex- 
trasolar planetary companions close to a bright parent star 
is a technically difficult task. With the development of so- 
phisticated image analysis and Adaptive Optics (AO) sys- 
tems on several modern large (8 m class) telescopes it is now 
possible. AO is required to both increase the peak intensity 
of the point-spread function (PSF) and to concentrate the 
photons which are scattered into a diffuse halo by the at- 
mosphere back into the diffraction limited core. Dedicated 
high -contrast imagi ng instruments, such as HiCAIO (Sub- 
aru, iHodapp etali 2008) . SPHERE (VLT , iBeuzit et all 
2008) and GPI (Gemini, Macintosh et all 2006), are de- 
signed to incorporate eXtreme AO (XAO) systems and so- 
phisticated coronagraphs to reject the light from the star 
whilst conserving the few photons from the angularly sepa- 
rated companion. Quasi-static speckles, mimicking the sig- 
nal of faint companions, are now li miting the detection ca- 
pabili ties of these ins truments (e.g. iFitzgerald fc Graham!. 
2006; Sou mmer et al.l . 2007). Quasi-static speckles in the fo- 
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cal plane are caused by non-common path errors and un- 
corrected aberrations in the primary mirror and other op- 
tical and mechanical components. If these aberrations were 
entirely deterministic they could be subtracted. However, 
quasi-static speckles are slowly varying aberrations making 
calibration difficult. 

The current most popular methods for quasi-static 
speckle reduction are PSF subtraction techniques such as, 
PSF estimation (e.g. Lafrcnicrc et al., 2007) angular dif- 
ferential imaging (ADI, M 

arois et all 2006 ). simultaneous 
spectral differential imaging (SSDI. ISmithl 1987) and, in 
the case of reflected li ght, polarimetric differential imag- 
ing (ISeager et al.l |2000T ). SSDI and PDI both require cer- 
tain properties from the target. ADI is more generic but 
if the speckles evolve during the observation the suppres- 
sion provided by the technique reduces dramatically. The 
temporal decorrelation timescale of these quasi-static speck- 
les is an important factor when estimating the performance 
of image subtraction techniques which have proved them- 
selves to be efficient a t Strehl ratios of the order of 20-40% 
l|Martinez et al.l l2012h . At high Strehl (>80%) the quasi- 
static speckles will become even more dominant as the PSF 
halo is reduced further. In this high Strehl regime the speckle 
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coherence timescale is unknown. iMartinez et al.l show that 
the quasi-static speckles become unstable over a timescale 
of a few seconds on a laboratory XAO test bench. It is 
thought that the evolution of the speckle pattern was pri- 
marily caused by temperature fluctuations and so, on a ded- 
icated instrument, this could be more controlled. 

Several other interesting and inventive techniques are 
also being developed in order to further enhance the prob- 
ability of detecting faint companions. iRibak fc Gladvsa 
l)2008l ) demonstrate it is possible to enhance the contrast by 
placing a rotating eccentric mask in the pupil plane. This 
breaks the symmetry between the telescope pupil and the 
focal plane causing the quasi-static speckles to move, the 
companion will however s t ay fix ed at the same location. 
iGladvsz fc Christoul ((2008, 2009) use the statistical distri- 
bution difference of on-axis and off-axis PSF to differentiate 
between real sources and speckles. It is likely that only the 
combination of several of these techniques in conjunction 
with AO and a coronagraph will result in the highest con- 
trast ratios in modern instrumentation. 

Here we propose to combine the idea of breaking the 
symmetry of the optical system a nd the sensor by c hang- 
ing the pupil function (similar to Ribak & Gladysz) with 
an adaptive pupil mask (APM) (|Osborn et al . 2009). The 
APM is positioned in the conjugate plane of the telescope 
pupil. This pupil mask consists of a number of independently 
controllable elements. The simplest design would be a seg- 
mented mirror where each segment can either reflect the 
light on-axis into the remaining optical system or off-axis 
into a baffle. This is also similar to the spec kle decorrela- 
tion (or phase-boiling) method (Saha 2002). The phase- 
boiling method involves adding additional phase aberra- 
tions to the optical path in order to force the quasi-static 
speckles to be more dynamic. This was later disproven by 
ISivaramakrishnan et al.l |2002f) as the original phase aberra- 
tions are still present and so the quasi-static speckles are also 
still present, simply hidden within a field of dynamic speck- 
les. Here we boil the speckles by introducing amplitude aber- 
rations instead of adding phase aberrations, an important 
difference as the interference pattern will no longer contain 
the same speckle pattern, but will actually be completely 
different. 

By changing the shape of the pupil we modify the 
diffraction pattern in the pupil plane. Light that once inter- 
fered constructively to form a quasi-static speckle at a given 
location will now not. However, speckles will be formed in 
other locations in the focal plane. The quasi-static speck- 
les will be forced to be dynamic, removing any dependance 
on speckle timescale. By changing the pupil function many 
times during an exposure these speckles will average out into 
a smooth PSF, albeit with an additional halo due to diffrac- 
tion through the pupil mask (Babinet's principle). This will 
mean that longer integration times are required for a com- 
panion to become visible with the same signal to noise ratio 
(SNR). However, the quasi-static speckles will be substan- 
tially reduced, reducing the complexity of the companion 
identification problem. 

In addition to the APM reducing the static speckles into 
a smooth halo, we can also expose after each configuration 
of the mask separately. As the quasi-static speckles are now 
dynamic we can use other image manipulation methods to 
further enhance the image. Dark speckle (DS) imaging, first 



proposed bv lLabevrid (|l995l ). is a technique designed for the 
detection of faint companions in the presence of dynamic 
speckles induced by the turbulent atmosphere. Briefly, after 
low order AO correction to focus the majority of the pho- 
tons in to the diffraction limited core and a coronagraph to 
reject the photons from the star there still remains some at- 
mospheric turbulence induced speckles in the PSF halo in 
short exposures. In this speckle halo, at some locations the 
wavefront will interfere destructively and result in a zero 
photon event. As the atmosphere evolves and traverses the 
telescope field of view, the position of these nulls will change 
in the focal plane. However, at the position of a faint com- 
panion the probability of a zero photon event is considerably 
lower. Therefore, by counting the number of times each pixel 
records a zero photon event in each short exposure we can 
generate a 'dark map' where the position of the companion 
will have a value lower than the rest of the image. Modern 
XAO systems are designed with target residual wavefront er- 
ror (WFE) of the order of a few nm and so there is actually 
very little in the way of residual atmospheric speckle. How- 
ever, as the quasi-static speckles are now dynamic we can 
use the dark speckle method to suppress the now dynamic 
quasi-static speckles. 

The DS method is intended to be used to enhance AO 
corrected images. Each DS exposure must be short otherwise 
the AO residual speckles will average, reducing the speckle 
nulls. Here we develop the technique for quasi-static speck- 
les, in which case the exposure time no longer needs to be 
short and we record the minimum value over a number of 
mask configurations. 

An additional advantage of the pupil mask is that it is 
a configurable device allowing it to act in several different 
observing modes. In 'adaptive' mode the mask can also be 
used to put a hard limit on the residual atmospheric WFE. 
Adaptive pupil masks like this have been shown to be able to 
reduce the PSF halo and actu ally increase PSF pe ak inten- 
sity despite removing photons l|Osborn et al.|[2009l ). making 
it useful in scenarios where only low-order AO or even no 
AO is available. The APM could be us ed in 'static ' mode as 
a non-redundant aperture mask (e.g. [Kopilovicb ], 1984) or 
partially-redundant aperture mask (e.g. iBuscher fc Hanifll . 
1993), used in many modern high contrast imaging instru- 
ments. If the APM is sufficiently high order it can also be 
used to emula te any binary shaped pu pil-pla ne mask coro- 
nagraph (e.g. iNisenson fc Papalioliosl . 2001; iKasdin et all 
2003). The configuration of these static masks can easily 
be changed to experiment with different positions, element 
sizes and configurations, a current area of active research. 

In section 2 we describe the simulation, in section 3 we 
introduce the APM and show how it will reduce the quasi- 
static speckles, section 4 describes application of the dark 
speckle method to APM, in section 5 we discuss the results 
and we conclude in section 6. 



2 SIMULATION 

2.1 Quasi-static speckles 

The combined mirror and optics abe r ration s are simulated 
using the method of ICavarroc et al. I (|2006l ) using the k 
(where n is the spatial frequency) power law in the spatial 
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Figure 1. The primary mirror phase aberrations and pupil in- 
cluding telescope secondary support spiders. The RMS error is 
one fifth of a wavelength. On the right is a log-scaled simu- 
lated image of a star with AO residuals, photon noise, read noise, 
sky background noise and quasi-static speckles. The lower panel 
shows the speckles. Plotted is the absolute difference between two 
images, one with quasi-static speckles and one without. 



power spectrum, as defined in iDuparree et all fH3). Fig- 
ure[T]shows a resultant static phase error, the corresponding 
image and an example of the speckles found by subtracting 
an image generated with quasi-static speckles from one with- 
out. 



2.2 Adaptive optics 

A Monte-Carlo simulation has been developed to test the 
concept. The simulation includes Poisson noise, sky back- 
ground noise (assuming 14 th magnitude) and 10e~ read 
noise. Residual adaptive optics phase aberrations are also 
included in the simulations. We assume a Strehl ratio of 
90% (a measure of AO residual error), consistent with pre- 
dictions for XAO systems. This residual wavefront error re- 
sults in a speckled intensity pattern around the diffraction 
limited core, over time this averages to a smooth halo. 

There are a few choices for simulating this AO corrected 
PSF: 

• We can assume a basic two part PSF, with a central 
diffraction limited core and a larger, in spatial extent, Gaus- 
sian halo with relative energy contributions corresponding 
to the desired Strehl. This would be the simplest case. How- 
ever, it is a very simple approximation and would result in 
a smooth estimate of the infinitely long exposure PSF. This 
will not include any low intensity and high spatial frequency 
residual speckles in the image, which are important to in- 
clude when developing a method to reduce such speckles. 

• We can use a fully analytical approach using statistical 
distributions to estimate the infinitely long exposure AO 
corrected PSF. As above, this will also result in a smooth 
PSF absent of residual speckles, but will be more accurate 
as it will include a deformable mirror model (subaperture 
sizes, mirror type etc.). 

• We can use a complete end-to-end Monte-Carlo simula- 



tion. This is more complicated and requires in depth mod- 
elling of the whole optical system (including wavefront sen- 
sors, reconstructor, mirror dynamics etc.) which are system 
dependant. This is a more complicated solution which is dif- 
ficult to develop, test and calibrate. 

• The last option is a combination of the two points 
above. We use a statistical expression to generate AO cor- 
rected phase screens - removing the need for wavefront sen- 
sors, reconstructor and deformable mirror dynamics, but re- 
taining the non-deterministic nature of the PSF. This will be 
an approximation which is not reflected by any real system 
but it will give an estimate of the AO corrected PSF con- 
taining residual PSF speckles averaged over whatever time 
scale we desire. 

Due to the instrument independent nature of this work we 
only need an approximation to an AO corrected PSF with 
the desired wavefront error. For this reason the corrected 
PSF is generated using the last method in the list above. 

The AO corrected PSF is estimated by summing over a 
number of independently realised instantaneous PSFs. Each 
of which is generated by a semi-analytical method of spa- 
tially filtering a von Karman phase screen with an AO trans- 
fer function and Fourier transforming to form an image. 

An AO system will reduce the spatial power spectrum at 
low spatial frequencies, the effe ct of which can b e modelled 
by a high pass filter, H(nd/2) (|Greenwoodlll978h . 



H(nd/2) 



2Ji M/2) V 
Kd/2 ) 



16 (2/nd) 2 Jl (nd/2) , (1) 



where d is the diameter of the subapertures and J„ is a 
Bessel function of the first kind of order n. The equation 
given is for a segmented mirror with tip/tilt and piston 
correction. This filter function only includes the deformable 
mirror fitting. Wavefront sensor errors (noise, non-linearity 
etc.), deformable mirror dynamics, latency and reconstruc- 
tor errors are not considered. 

The AO residual phase spectral density, <3?ao is then 
given by, 



*A0(k) = #(rcd/2)$ a t m o S («) 



(2) 



where $ a t mos (K) is the von Karman spatial phase power 
spectrum ($atmo S (t) = r^ 5/:i (k 2 + 1/Lq) 11/6 , r is the 
Fried parameter - a measure of the strength of the turbulence 
and I/o is the outer scale of the atmospheric turbulence) . 

The AO co rrected phase screen is then given by 
l|Ellerbroekll2002l ), 



0.1517 



]V ) ' StF[v^(^)(x(«) + «(«))] 



V2 \ro J 

(3) 

where W is the width of the phase screen and (x( K ) + *x'( K )) 
is a randomly generated repeatable white noise field. 

The filtered optical turbulence phase, 0ao , is added to 
the static phase aberrations, Static- This assumes that the 
two are independent. In fact, the AO system will attempt 
to correct some of the static aberrations and so ^static is 
also spatially filtered by the AO filter function. Therefore, 
$ AO, static = jy"(Kd/2)$ sta tic, can be used to replace <E>ao(k) 
in equation[3]to generate the quasi static phase aberrations. 



4 J. Osborn 



The PSF is given by, 

N 

Hp, f) = £ l-F LP(0 ex P (-i(^tatic(f, »?) + 0ao(C, t)] I 2 , 



(4) 

where T is the fourier transform operator, P(£) is the pupil 
function and is equal to 1 when £ < D/2 and otherwise, 
D is the diameter of the telescope aperture, p and 9 are the 
polar co-ordinates in the focal plane, £ and 77 are the pupil 
plane variables in polar co-ordinate space, iV is the number 
of simulation iterations and t denotes the time variable. 




Figure 2. Two individual consecutive images after APM (left and 
centre). The speckle pattern in each image is different as shown 
by the modulus of the difference of the images (i.e ii — Ta|, where 
Ii and I2 are the intensity patterns of the two frames) in the right 
panel. In this example we block a random 15% of the pupil. 



2.3 Simulation parameters 

The process is repeated and a long exposure image is built up 
for as many iterations as required. The telescope diameter 
is 8 m, d = 0.18 m, ro = 0.12 m, Ln = 30 m, the observing 
wavelength is 1.6 /im with a bandwidth of 0.23 /im (H-band), 
the exposure time is 30 s and the pixel scale is 0.25A/D. 

No coronagraph is included in the simulations, this is 
because modern coronagraphic techniques are as numerous 
as they are complicated. This means that any achieved con- 
trast ratios can not be directly compared to those from any 
high contrast imaging instruments. Instead we compare to 
a standard image with AO but without any other form of 
image manipulation. 

Figure [T] (top right) shows the simulated focal plane 
image of a bright target and nine fainter companions. The 
companions are distributed from the centre to the bottom 
of the image. The central star has a magnitude of 4 and 
the companions have magnitudes 14- 22, reducing by one 
magnitude for each step (2.5A/D) outwards. 



3 ADAPTIVE PUPIL MASKING 

A re-configurable mask is used to block the light of a cho- 
sen fraction, /, of the pupil. This is done with a spatial 
light modulator with 44 x 44 elements. In each iteration a 
random selection of these elements are flipped to block an 
area of the pupil. The pupil function therefore changes in 
each iteration. As the image can be thought of as the in- 
terference pattern of the pupil function multiplied by the 
complex amplitude of any phase aberrations, the image in 
the focal plane will therefore also be modified. The random 
nature of the APM is important as this means that it is com- 
pletely independent from any other optical component, it is 
modular. No additional information or optical components 
(e.g. wavefront sensors) are required. 
The image is now given by, 

N 

i( P ,e) = Y l \J r [ M (^,t)P(0 


exp (-i(0stattc(£, rj) + 0ao(£, »?,<)))] | 2 , (5) 

where M(£, r/, t) is the time dependant mask function. In 
the simulation presented here, we block a random 15% of 
the telescope aperture in each iteration. The actual fraction 
can be chosen and optimised by the user depending on sys- 
tem specific parameters. Here we choose 15% as this is con- 
sistent with the through put of many coronagraphic systems 
l|Ribak fc Gladvszll2008r i. 
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Figure 3. The sum of 1000 images. In this example we block a 
random 15% of the pupil. 



In the modified images the quasi-static speckles are no 
longer quasi-static but are actually very different in every 
iteration. Figure [5] shows two consecutive focal images and 
the modulus of the difference. In each individual frame the 
speckle pattern is completely different. No atmospheric tur- 
bulence was included to generate these plots and so the 
speckles movement is entirely due to the changing pupil 
mask. Figure |3 shows the sum of 1000 frames. We can define 
a metric for the magnitude of the quasi-static speckle as the 
normalised azimuthal variance, 



\I(p,0)-{I(p))\ 
(HP)) 



-de, 



(6) 



where (I(p)) denotes the expected intensity of the PSF at 
radius p. The metric must be normalised by the azimuthal 
average or the variance will be biased by the radial inten- 
sity of the image. Figure [4] shows the azimuthal variance for 
the original image and the APM image. We see that the 
radial intensity variance for the original image is approxi- 
mately constant for all field angles, the speckle is dominant 
over the diffraction halo. The azimuthal variance for the 
APM image continues to reduce with separation. At small 
angles the diffraction pattern from the central star is in- 
tense and the speckles are 'pinned' to the first diffraction 
rings (|Bloemhof et al.ll200ll ). even with the APM. At larger 
angles the speckles are more free to move and so we achieve 
a greater suppression. 

When we sum the images over a number of iterations, 
each with a different pupil pattern, the quasi-static speck- 
les average out to a smooth PSF. However, as we are using 
square blocking elements Babinet's principle dictates that 
we can expect the diffraction limited PSF to be a super- 
position of a square diffraction pattern of the blocking el- 
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Figure 4. Azimuthal variance as a function of radius for the 
original image (black line) and the APM image (dashed line). 



ements and the circular diffraction pattern from the pupil. 
The diffraction limited PSF in polar co-ordinates is therefore 
given by, 



PSF(p, 



1 



(1 



2J 1 (pL>/2) 2 2J 1 (apD/2)\ 



pD/2 



(1 - /) + 2a sine 



pa sin(0) 



apD/2 J 
pacas(e)\ j. ^ 



where p is related to the pupil plane parameters by p — 
27r£/A/i, A is the wavelength of the light and fi is the 
effective focal length of the optical system, a is the fractional 
radius of the central obscuration and a is the length of one 
side of one of the masking elements. 

As the mask elements are smaller than the pupil, the 
diffraction pattern is broader, increasing the halo around 
the PSF. Smaller blocking elements will result in a broader 
diffraction halo and a greater blocking fraction will result in 
a diffraction halo with a greater fraction of the total energy. 

Although the sum of the individual frames does show 
a reduction in the quasi-static speckles, the disadvantage 
would be the decreased the signal to noise ratio (SNR) due 
to the increased halo intensity and reduced peak intensity. 
This means that we would have to integrate for longer to 
achieve the same SNR, although the quasi-static speckles 
will no longer form a fundamental limit - we are now photon 
noise limited. 

The theoretical form of the APM PSF is derived in 
orn et all l|2010r i. Here we briefly review the theoretical 
structure of the PSF for the case of random pupil blocking. 
The PSF, assuming on-axis observations, can be estimated 
from the modulation transfer function of the residual atmo- 
spheric phase aberrations and of the telescope aperture. The 
PSF is given by, 



PSF = [F[MTF ao x MTF tcl ]) x (1 - /), 



(8) 



where MTF ao is the atmospheric modulation transfer func- 
tion and MTF tc \ is the telescope modulation transfer func- 
tion. As the MTF is defined for unit total intensity the PSF 
must be scaled by (1 — /). We can see that the total intensity 
of the image is reduced by / but the change in peak inten- 
sity is not so obvious. It has been shown that with careful 
pupil blocking the p eak intensity can actually be increased 



(Osb orn et al 



_e_P' 

120091 ). However, here we choose to randomly 



block the pupil. In this case we would expect the peak in- 
tensity to be reduced. Peak intensity is equal to the integral 



of the MTF, 
Io = 



{MTF ao x MTF tcl )du{l - /), 



(9) 



where v is the spatial frequency in the focal plane and is 
related to the separation in the pupil plane, r by r = X/lV- 
The MTF is related to the phase structure function D by, 



MTF ao{v) = exp(-0.5D A o(f)), 
and (|Rao et al.ll200oh 



D AO {v) 



47T / [I — Jo(Kv)\$AOKdK. 
'0 



(10) 



(11) 



Therefore, 
MTF ao iy) 



exp 



^ — 0.5 X 47T J [1 — Jo(K^)]$AOKdK 



(12) 



MTFtd is given by the autocorrelation function of the aper- 
ture function. This will be the product of the pupil function 
and the mask, 



and so, 



MTFtd(v) =F[\W(t,r},t)\ 2 ] 



(13) 



(14) 



The mask will act to reduce MTF tc i due to the lower fill 
factor in the pupil (less redundency in Fourier baselines), 
reducing the peak intensity of the image. Using the above 
we can write equation [H] as, 



PSF = T 



exp ( —0.5 x 4-7T / [1 — Jo^nu^^AOKdn 



X?[\W(t,Tl,t)\ 



*(!-/), (15) 



and equation [9] as, 
I = I I exp 



0.5 x 4n [1 — Jo(Kv)]&AOKdn 



xP[|VK(£,r?,t)| 2 ])aVx(l-/), (16) 

This is solved numerically and for a random blocking of 15% 
we find that the peak intensity is reduced by ~27%. To take 
a specific example, if we would like to observe a companion 
at 5X/D the normalised halo would be increased by ~16% 
in comparison to the non-blocked image. 

The exposure time, to, to achieve a given SNR assuming 
we are limited by the photon noise of the halo is given by, 



to = SNR 2 c/SAe, 



(17) 



where S is the signal from the planet, A the telescope area, 
e is the efficiency of the optical system and c is the contrast 
ratio between the intensity at the position of the image and 
the companion signal, c = n 7 ,h(p = 5A/D)/n 7 , comp (p = 0), 
where n 7i h(p = 5X/D) is the number of photons in the halo 
at the position, p = 5X/D and n 7 , comp (p = 0) is the peak 
number of photons from the companion. 

As the mask will increase the PSF halo and reduce the 
companion signal the contrast will be increased, requiring 
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Figure 5. Dark speckle image. Plotted is the minimum count of 
each pixel over every iteration of the mask. Eight of the nine com- 
panions are now clearly visible. The ninth, the closest, companion 
is still hidden in the halo of the bright star. 



longer exposure times to conserve the SNR. For the ratio 
to/tAPM this simplifies to, 



to/tAPM = (1 — f)c/CAPM, 



(18) 



where tAPM and capm is the exposure time and contrast 
ratio of the APM image respectively. Using the figures found 
numerically above we see that the contrast is increased by 
a factor of 1.59 for a companion at 5X/D. The exposure 
time to achieve the same SNR in the noiseless case would 
be 1.35 times longer than the non-blocked image. This is 
confirmed by the Monte-Carlo simulation with a value of 
1.37. With noise (background, read and shot noise) the value 
is increased to 1.47. 



4 DARK SPECKLE IMAGING 

As we have forced the quasi-static speckle pattern to be 
dynamic we can use the individual frames to increase the 
image quality. One suggestion, developed here, would be to 
use the dark speckle analysis of lLabevriel fl995). 

In each frame the speckle pattern will have minima in 
intensity. The location of these nulls will vary as the pupil 
mask changes shape. In the position of the faint companion 
the probability of seeing a speckle null is considerably lower 
than anywhere else. 

The proposed technique here is slightly different to that 
of Labeyrie. Dark speckle imaging is not applicable to long 
exposures as the atmospheric speckles will average over time, 
reducing the depth of the nulls. However, here we are not 
necessarily restricted to very short exposures and hence the 
read out noise, sky background and AO residuals will mean 
that the probability of measuring zero photons over an arbi- 
trary time period is low. Instead we record the lowest value 
of each pixel over all of the iterations of mask configuration. 
We have chosen to expose 1000 iterations at 3 s exposure 
time each, resulting in a total observing time of approxi- 
mately 50 minutes. Figure [5] shows the dark speckle image. 
We see that eight of the companions are now easily visible. 
The companion closest to the bright star (at 2.5X/D) is still 
hidden in the halo from the bright star. 

A disadvantage of this technique is that the intensity 
at the location of the companion will be the minimum flux 
measured from the companion over an exposure time of one 
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Figure 6. Focal images assuming a perfect coronagraph. The 
images shown are for the original image (left) and the DS image 
(right). 
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(a) 
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(b) 
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(c) 

Figure 7. Focal images including nine faint companions aligned 
from image centre to the bottom. The images shown are for the 
case with no manipulation (a), long exposure APM image (b) and 
APM + DS image (c). All images are now normalised relative to 
the peak intensity of the original image. 



iteration. This will have the expected value of the intensity 
of the companion minus the shot noise. For this technique 
to work a large number of readouts are required, although 
it should be noted that as we are only recording the min- 
imum value of each pixel we are essentially only including 
one readout per pixel. 

A coronagraph would reject a large fraction of the pho- 
tons from the parent star and improve the results at small 
angular separations. For example, if we assume a perfect 
coronagraph (i.e. one which removes the diffraction limited 
PSF but not the quasi-static speckles, noise or any diffrac- 
tion effects due to the mask) then figure [6] shows the original 
image and the DS image. All companions are now visible 
with very little speckle. 



5 DISCUSSION 

For comparison, figure [7] shows the final images including 
the nine faint companions for the case of no manipulation, 
long exposure APM and APM with DS image analysis. We 
can examine radial cuts of the final images to compare per- 
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Figure 8. Radial cut along the direction on the simulated com- 
panions (black) and azimuthal average (grey) for the thresholded 
APM+DS image, long exposure AP image and the original image 
(from top to bottom). The vertical lines indicate the positions of 
the companions in the focal plane. 
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Figure 9. Radial cut along the direction on the simulated com- 
panions (black) and azimuthal average (grey) for the APM+DS 
image (dashed line), the long exposure APM image (dotted line) 
and the original image (thin line). The vertical lines indicate the 
positions of the companions in the focal plane. 

formance. In figure [8] we show a radial cut in the direction 
of the companions and the azimuthal average for the final 
APM+DS image (top), long exposure APM image (middle) 
and the original image (bottom). The vertical lines indicate 
the position of the faint companion. Figure [9] shows the ra- 
dial profile of the original image the long exposure APM 
image and the APM+DS image on the same axis. 

We can see that the effect of the companions is visi- 
ble in the original and long exposure images, but only the 
APM+DS method has been able to isolate all but the closest 
companion successfully in this case. At separations greater 
than approximately 5X/D the contrast between the PSF 
peak and the halo is increased by the DS method by a factor 
of approximately 100, corresponding to 5 magnitudes. It is 



clear that not only have the speckles been suppressed but 
also the bright star diffraction pattern is also suppressed. 
This is because the speckles nulls are lower than the average 
intensity of the PSF halo for any given location. We are now 
able to detect companions 11 magnitudes fainter than the 
star at separations of 5X/D and 18 magnitudes fainter at 
22.5A/D. The inner working angle would be improved with 
the use of a coronagraph or other PSF subtraction methods 
which is not included here. 

The performance of the technique is dependant on the 
parameters selected and on the properties of the system we 
wish to observe. Detection is modified if we increase the frac- 
tion that the APM blocks, the number of iterations or the 
Strehl ratio of the AO system, of which, only the blocked 
fraction of the mask and the number of iterations are con- 
trolled by the observer. If we increase the fraction which is 
blocked then we generate more speckle in the focal plane, 
which means quicker averaging and more nulls. However, it 
will also mean reducing the throughput of the system, there- 
fore this will depend on the system to be observed. In the 
simulation we blocked 15% of the pupil. The number of iter- 
ations of the mask is also important. The more iterations the 
greater number of configurations of the mask are used. This 
means more averaging in the long exposure APM image and 
means more nulls are measured in the DS image, increasing 
detection. The number of iterations is a balance between the 
amount of time available for an observation and the expo- 
sure time of each iteration. In order to receive a significant 
number of photons from the companion it may be necessary 
to increase the exposure time, this will limit the possible 
number of iterations. In the simulations we assumed 1000 
iteration of 3 seconds each. It should be noted that there is 
no limit to the exposure time of each iteration. 



6 CONCLUSIONS 

Atmosphere induced phase aberrations (with or without 
AO) causes a halo of speckles to form around the PSF. This 
source of noise will average over time and, given enough 
time, faint companions can be observed above this halo. Cur- 
rent high contrast imaging is limited by quasi-static speck- 
les caused by the optics and structure of the telescope. This 
is because these quasi-static speckles do not average. They 
appear in the image as potential false-positive candidates 
for faint companions and are difficult to distinguish. Here 
we have presented a technique to turn these static speck- 
les into dynamic speckles. This means that over time these 
static speckles will also average into a broad halo. The shape 
and magnitude of this halo will depend upon the geometry 
of each mask element and the fraction of the pupil that is 
blocked. This means that over time the PSF will converge to 
a smooth form allowing the companions to be seen above. 
However, as the APM adds more energy to the halo and 
the throughput of the system is reduced by a fraction equal 
to the fraction of the pupil which is blocked, the SNR is 
reduced. We would need to observe for a longer period of 
time to collect enough photons from the companion to be 
seen above the halo. For an example case of a companion 
at 5A/.D we would need to observe for 1.35 times longer 
(1.47 defined by the simulation including noise) to achieve 
the same SNR. In this case the we do not need to expose the 
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CCD between each mask state and so the mask frequency 
can be high allowing many mask states for an arbitrarily 
long exposure. 

In addition to the smoother PSF we can also use other 
conventional image manipulation techniques which are nor- 
mally used for the dynamic speckling caused by atmospheric 
turbulence. Here we show the effect of implementing an 
adaptation of Labeyrie's dark speckle analysis. This is a sim- 
ple technique where we record the minimum value of each 
pixel over all of the iterations of the mask. 

Using this technique in a Monte-Carlo simulation we 
can detect faint companions with a higher magnitude dif- 
ference from the central star. We find that at separations 
greater than a few X/D the PSF halo count is reduced by 
a factor of approximately 100, corresponding to a contrast 
increase of 5 magnitudes. 

Due to the APM the diffraction rings from the bright 
star also move, this means that the diffraction pattern from 
the bright star is also suppressed by the APM when used 
with the DS imaging analysis. At small inner working an- 
gles there is still some confusion between the companions 
and that of the parent star, due to the pinned speckles in the 
bright diffraction rings. However, we are now able to detect 
companions 11 magnitudes fainter than the star at separa- 
tions of 5X/D and up to 18 magnitudes fainter at 22.5A/D. 
The inner working angle would be improved with the use of 
a coronagraph or other PSF subtraction methods which is 
not included here. In order to perform the DS analysis short 
exposure images will have to be recorded. The length of the 
short exposures is arbitrary but the optimal will depend on 
the number of iterations required during the exposure and 
the magnitude of the target. 

The noise attenuation found here is comparable with 
the a ttenuation from the H igh-Order Testbench (HOT) with 
ADI (|Martinez et alj[2012f ). The advantage of this technique 
over ADI is that it does not rely on PSF subtraction and 
is therefore insensitive to any changes in the quasi-static 
pattern. It has the advantage over PDI and SSDI that it 
also insensitive to the companion properties, i.e. it does not 
require the companion to have a specific emission spectrum 
or the light to be reflected and hence partially polarised. 
The disadvantage is the reduced throughput. 

APM would benefit from collaboration with other tech- 
niques. Currently the inner working angle is limited by the 
diffraction pattern of the primary star. A coronagraph would 
reject light from the central star resulting in even higher 
achievable contrast ratios. Also, no additional image ma- 
nipulation was performed. No attempt at PSF subtraction 
was made. The dynamic speckles will average into the pre- 
dictable smooth long exposure PSF, which could be sub- 
tracted. Combination with other techniques would improve 
the noise attenuation and reduce the probability of false- 
positives. 

It is also important to note that by including a high 
order APM, in addition to quasi-static speckle removal, it 
would also be possible to reduce the residual WFE, emulate 
a non-redundant (or partially redundant) aperture mask and 
a binary shaped pupil plane coronagraph, all of which would 
be completely and easily re-configurable. 
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